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The relativistic jet of Cygnus X-3 in gamma rays 
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ABSTRACT 

High energy gamma-rays have been detected from Cyg X-3, a system composed of a 
Wolf-Rayet star and a black hole or neutron star. The gamma-ray emission is linked 
to the radio emission from the jet launched in the system. The flux is modulated with 
the 4.8 hr orbital period, as expected if high energy electrons are upscattering photons 
emitted by the Wolf-Rayet star to gamma-ray energies. This modulation is computed 
assuming that high energy electrons are located at some distance along a relativistic 
jet of arbitrary orientation. Modelling shows that the jet must be inclined and that the 
gamma ray emitting electrons cannot be located within the system. This is consistent 
with the idea that the electrons gain energy where the jet is recollimated by the stellar 
wind pressure and forms a shock. Jet precession should strongly affect the gamma-ray 
modulation shape at different epochs. The power in non-thermal electrons represents 
a small fraction of the Eddington luminosity only if the inclination is low i.e. if the 
compact object is a black hole. 

Key words: radiation mechanisms: non-thermal — stars: individual (Cygnus X-3) 
— ISM: jets and outflows — gamma rays: theory — X-rays: binaries 



1 INTRODUCTION 

Cyg X-3 is a high-mass X-ray binary composed of a com- 
pact object in a 4.8 hr orbit arou nd a Wolf-Rayet (WR) star 



at a distance of about 7 kpc (se e Bonnet-Bidaud fc ChardTnl 
! Il988l : Ivan Kerkwiik et al.Nl996l ; iLing et alj|2009t and refer- 



ences therein). The system is a bright X-ray source with 
Lx ~ 10 38 erg s _1 . Cyg X-3 is also well-known for radio flar- 
ing (up to 20 Jy) w hen the source has a soft X-ray spectra 
i|Szostek et al.ll2008r ). The radio source is resolved into a rel- 
ativistic jet with an expansion speed of 0.3-0. 7c. The strong 
stellar wind from the WR companion (M w « 10 _5 Af s yr- 1 , 
v w ~ 1000 km s _1 ) has a major impact on the environ- 
ment of the high-energy source. Scattering in the wind is 
probably responsible for washing out rapid X-ray variabil- 
ity timescales and also for modulating the X-ray emission. It 
acts as a veil that has made it difficult to identify the nature 
of the compact object, black hole or neutron star. Despite 
the differences caused by the WR wind, Cyg X-3 is firmly 
established as a trademark accreting binary with relativistic 
jet i.e. a microquasar. 

The AGILE and the Fermi Gamma-ray Space Tele- 
scope collaborations have recently reported the detection of 
high-energy gamm a rays (HE, >100 MeV) from Cyg X-3 
i|Tavani et al.ll2009l : IXbdo et alj|2009l h The identification is 
firm because the detections occur exclusively when Cyg X-3 
is flaring in radio and because Fermi observations show the 
HE gamma-ray flux is modulated with the orbital period. 
The gamma-ray modulation is almost in anti-phase with 



the X-ray modulation, with the gamma-ray minimum oc- 
curring about 0.3-0.4 in phase after X-ray minimum. The 
modulation amplitude is close to 100% after background 
subtraction. The spectrum is consistent with a power law 
F v ~ v~ a with a = 1.7. The luminosity above 100 MeV is 
a few 10 3B (d/7 kpc) 2 erg s _1 . 

Inverse Compton (IC) scattering of photons from the 
WR star on high energy electrons is a natural candidate to 
explain the gamma-ray emission. The high temperature of 
the WR star {R* fti 1 Rq, T* « 10 5 K) and tight orbit 
(d ~ 3 10 11 cm) imply that the radiation density in pho- 
tons from the star is u* « 10 s erg cm -3 at the location of 
the compact object, which is at least an order-of-magnitudc 
higher than any other X-ray binary. Electrons with Lorentz 
factors of a few 10 3 upscatter 20 eV stellar photons above 
100 MeV very efficiently in such a radiation field. IC scat- 
tering directly produces a modulation of the flux because of 
the orbital motion. The maximum occurs when stellar pho- 
tons are backscattered towards the observer. The accretion 
disc can also provide seed photons if the HE electrons are 
close enough. This does not lead to a modulation unless the 
HE electrons - disk geometry seen by th e observer changes 
with orbital phase (Meszaros et al.lll977l ). Pion production 
is possible if there are high energy protons. However, even 
in this dense environment, it is less efficient than IC so that 
its energy requirements are higher. 

The link between gamma-ray and radio flares suggests 
that the HE electrons are located in the relativistic jet. 
Observations of knots in active galactic nuclei show that 
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Figure 1. Geometry of the jet model. The scattering electrons 
are situated at symmetric locations in a jet with relativistic speed 
0. The seed photon source is the star. 



particles may be accelerated at specifi c locations along th e 
jet, linked e.g. to recollimation shocks (jStawarz et al.| [2006). 
Assuming the electrons mainly upscatter stellar photons at 
some location along the jet, the expected IC emission will 
depend upon the distance to the star, the bulk velocity of 
the jet and its orientation. This orientation is not necessarily 
perpendicular to the orbital plane if e.g. the inner accretion 
disc is warped or it depends on the black hole spin axis. 
However, the jet orientation is fixed as seen by the observer 
(changing only if the jet precesses). 

The goal here is to test quantitatively whether the 
Fermi gamma-ray modulation can be reproduced in this 
framework and to see if constraints can be derived on the 
jet parameters. 



2 JET INVERSE COMPTON EMISSION 
2.1 Emission spectrum 

The HE electrons are assumed to be located at a distance 
H from the compact object along a jet with a bulk velocity 
/3 — v/c (Fig. [TJ. The stellar emission is approximated as 
a point-like blackbody of temperature T* and luminosity 
AhR 2 <jsbT+ . The electron Lorentz factors 7 e are distributed 
as a power-law dN e = K e n f^r p d n f e . In the Thompson regime, 
the inverse Compton e mission spectrum a t a photon energy 
e (in ergs) is given by (|Dubus et alj|201oh 



(1 - e*.e obs ) c 



x ^obs 



(1) 



where: the flux index is related to the electron power law 
index through a — (p — l)/2, R is the distance from the 
star to the electron location; and e Q b s are unit vectors 
along, respectively, the star-to-electrons and the electrons- 
to-observer directions; 



(1-/3 2 ) 



2\l/2 



(1 — /3e bs-ejet) 



(2) 



defined the Doppler boost of the jet, ej ct being the unit 
vector along the jet direction; C(p) is given by 



C(p) = 



h 3 c 3 



(p 2 + 4 P + n)r(^)c(^) 

(p+l)(p + 3)(p + 5) 



(3) 



with F the gamma function and £ the Riemann function. 
This formula is valid in the Thompson regime, that is when 
7 e eo < m e c 2 where eo is the characteristic energy of the seed 
photons. For a blackbody with T* = 10 5 as in Cyg X-3, eo ~ 
2.7kT+ « 23 eV so the limit occurs for 7 e « 2 10 4 (neglecting 
the Doppler boost). IC emission from 100 MeV to a few GeV 
(the relevant Fermi range) occurs in the Thompson regime. 

The model geometry is shown in Fig. [T] The jet has 
an azimuth 9j and polar angle cj>j (=0 when perpendicular 
to orbital plane). With the origin set at the location of the 
WR star, 



R 2 =d 2 +H 2 + 2dff(e c .eje 



(4) 



where e c is the unit vector along the star to compact object 
direction, and the unit vectors are given by 



e 4 
e je t 



(de c + He jet )/R 



= (cos i 



sm ( 



sin U\ sm < 



e c = (cos 9, sin 9, 0) 
e bs = (0, — sini, cosi) 



(5) 



with 9 the true anomaly, d the orbital separation and i the 
inclination. Here, the true anomaly is defined so that 9 — 
±7r/2 at conjunctions. 

2.2 Main properties 

The inverse Compton emission has an orbital modulation 
because of the dependence of e c on the true anomaly (= 
orbital phase for a circular orbit). Developing dFic/d9 — 0, 
the emission maximum and minimum along the orbit verify: 

(a+l)(e c xe obs ).e z = — ((a + 3)e*.e obs - 2) (e c xe jot ).e z (6) 
R 

If H <C d, or if the jet is perpendicular to the orbital plane, 
then the maxima and minima are at conjunctions as outlined 
in §1. Otherwise, they occur at orbital phases that can be 
very different. 

The IC flux will be equal to zero if e*.e b s = 1 some- 
where along the orbit. Having a 100% modulation can be 
translated into a necessary condition on H for given i, d, 
<j)j and 9j. Similarly, although the seed photon density de- 
creases with H , the maximum of the IC flux for a given jet 
geometry does not necessarily occur for H=0 because of the 
dependence of e* on H. 

The jet speed only appears in Dobs and e b a .ej e t is con- 
stant along the orbit: changing (3 will only impact the flux 
normalisation and not the shape of the modulation. The 
maximum flux occurs when j3 = e oba .ej e t. Emission from a 
jet oriented away from the observer will always be weak for 
highly relativistic speeds because of the deboost. 



3 APPLICATION TO CYG X-3 

The observed modulation is plotted in Figure The back- 
ground level in diffuse gamma rays of 3.6 10 ~ 6 ph cm" 



s 1 was subtracted to the Fermi lightcurve (|Abdo et al.l 
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requires that 7 e ,mm ^ 1000 so P e is a lower limit on the 
non-thermal power. 

Good fits can be obtained for both Ol (Xmm = 2.7 
for 12 data points - 5 parameters = 7 degrees of freedom) 
and 02 (xLn = 4.2). The best model for Ol is plotted in 
Figured It has j3 = 0.41, H = 8 10 u cm, fa = 39°, $j = 
319°, P e = 10 38 erg s" 1 . The 90% confidence range for the 
par ameters was determin ed by adding 9.24 to the minimum 
\ 2 l|Lampton et all 1 19761 ) . Only models that had P e lower 
than the Eddington luminosity L Edd « 10 38 (M 1 /M Q ) erg 
s _1 were kept. Besides being physically implausible, models 
with larger P e are associated with high values of /3 or large 
H. The high P e then compensates for Doppler deboosting 
or low IC efficiency (see §3.3). 



Figure 2. Model fits to the observed > 100 MeV gamma-ray 
modulation in Cyg X-3. Conjunctions are at phases 0.25 and 0.75 
for the conventions adopted in this work. The models shown as- 
sumed an orbit with a black hole (Ol). The best model is shown 
with a black solid line. A model with = is shown with a grey 
solid line. The model with minimum P e (3 10 33 erg s _1 ) is shown 
with a grey dashed lines. All of these models are statistically ac- 
ceptable fits to the data (see §3 for details). 



2009). There is not absolute phasing of the orbit of Cyg X- 
3. The Fermi observations have been phased so that the 
well-defined minimum X-ray flux occurs at superior con- 
junction i.e. phase 0.25 with the conventions adopted in 
this paper (Fig. [T]). This phasing is justified if the X-ray 
modu lation is due to Thompson scattering in the stellar 
wind l|Pringlel Il974l ). It is independently supported by in- 
frared spectroscopy (|Hanson et alJl200Ch . 

The orbital parameters of Cyg X-3 are not determined 
precisely (|Hanson et al.ll2000l ; IVilhu et all [ 2009) so two ex- 
treme solutions are adopted following ISzostek fc Zdziarskil 
(120081 1. Orbit 1 (Ol) has a Mi=20 M Q black hole around a 
50 Mq WR star of radius 2.3 Rq and is seen with an incli- 
nation of 30°. Orbit 2 (02) has a Mi=1.4 Mq neutron star 
around a 5 Mq WR star of radius 0.6 Rq with i — 70°. The 
Fermi spectrum a — 1.7 sets the electron power-law index 
p = 4.4. The emission arise from two symmetric sites: the 
jet and the counterjet. The counterjet has <f> c j = tt + fa. 



3.1 Parameter exploration 

The jet is parametrised by /3, H, fa, 9j and K e . The expected 
modulation in the Fermi band is calculated using the equa- 
tion in §2 for the jet and the counterjet. The evaluation of 
Eq.[T]is very fast and allows an exhaustive exploration of the 
parameter space. The jet angle fa was varied between and 
7r/2 ; 8j varied between and 2ir. The emission height H 
was varied between O.Old and lOOd in logarithmic steps (d is 
the orbital separation). The jet speed /3 was varied linearly 
from to 0.99 (bulk Lorentz factor w 7). 

The model K e is adjusted to minimize the x 2 goodness- 
of-fit to the observed modulation. The normalisation K e is 
converted into a power in HE electrons P e assuming a dis- 
tance of 7 kpc and a minimum HE electron Lorentz factor 
7c, min = 1000. P e is highly sensitive to 7 e , m in because of the 
very steep electron spectrum. IC emission above 100 MeV 



3.2 Jet orientation 

Figure |3] shows the distributions of /3, H, fa and 8j for the 
black hole case (Ol). The figure also shows the distributions 
for various limits on P e . In all cases, the HE electrons dis- 
tance H is between 0.5 and 30 times the orbital separation 
(i.e. between 2 10 11 and 10 13 cm). A location very close to 
the compact object is excluded. The orientation of the jet 
is constrained to be 20° < fa < 80° with a preference for 
values comparable to the system inclination (i — 30°). A jet 
perpendicular to the orbital plane does not fit the data. The 
azimuth 6j is less constrained: there is a well defined peak 
in the distribution (bottom panel, Fig. [3} but, contrary to 
H or fa , there are good models all over the range even if in 
small numbers (not visible on a linear scale). 

Moderate relativistic speeds /3 are favoured but this is 
not strongly constrained. The speed is closely linked to the 
power in HE electrons. There is a tendency to have lower 
values of /? when the allowed P e gets smaller, accompanied 
by a smaller H. A model in the 90% confidence region with 
/3=0 is shown in Figure [2] It has X 2 = 7.1, H = 7 10 n cm, 
fa = 31°, 6j = 9°, P e = 2 10 37 erg sT 1 . This trend on /3 
reverses for low values of P e < 0.001 I/Edd- These do not 
appear in Figure [3] as there are comparatively very few such 
models. The minimum P e in the 90% confidence region is 4 
10 33 erg s _1 , a very modest fraction of LeoM- This model is 
also shown in Figure [2] It has x 2 = H-3, /3 = 0.99, H = 



have fa w i and 9j s=s= 



90°: they are almost aligned with 
the observer (ej e t.e b s ~ 1) at superior conjunction. The 
slight difference in 8j accounts for the phase difference of the 
maximum. Here, Doppler boosting compensates for the low 
P e . There is some degeneracy between the two parameters 
up to some (large) value of the Lorentz factor ~ 20 where 
good models cannot be found anymore. These are effectively 
microblazar models. 

The constraints in the neutron star case (orbit 02, not 
shown here) are similar. The jet orientation is well con- 



strained with 25° < 



< 65°, 



-60° < 0j < 



-10° and 



2 10 11 cm < H < 6 10 11 cm (H/d from 1 to 3), comparable 
to the values found with Ol. However, in all cases j3 is < 
0.2. Interestingly, P e is constrained to be rather large with 
P e > 0.2I/Edd (about 3 10 37 erg s _1 ). The large inclination 
(70°) required for a neutron star primary is the reason for 
the difference with the black hole case. Arbitrarily setting 
i — 30° with the orbit 02 gives results for /3 and P e that are 
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Figure 3. Distribution of jet parameters for models in the 90% 
confidence region given by \ 2 statistics. Orbit Ol (20 Mq black 
hole, i=30°) is assumed. The various regions correspond to a 
power in high energy electrons P e s£ -L E dd (light grey), ^ O.lLEdd 
(grey), < 0.01L E dd (dark grey). Here, L Edd is 2 10 39 erg s _1 . 



consistent with those of Ol. Large inclinations do not allow 
good fits for small values of P e or large values of /3. 

These results were obtained for a steep power-law dis- 
tribution of electrons with an index p = 4.4, because of the 
soft gamma-ray flux index and the assumption of Thomp- 
son scattering. Taking p = 2 or p — 3 does not affect the 
conclusions. A few tests cal culations using the full IC cross 
section (done as explained in lDubus et al.ll2010l ) showed that 
a slightly harder electron index (p « 4) is required to match 
the spectrum. Again, this does not change the results. The 
steep spectrum may not directly reflect an electron power- 
law distribution but represent the best fit to e.g. a cutoff 
in the 100 MeV - 1 GeV range. To test this, a lightcurve 
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Figure 4. Impact of jet precession on the gamma-ray lightcurve 
for the best-fit model shown in Figure [2] The jet azimuth 0j is 
rotated in steps of 72° from its best fit value of 319°, with lighter 
lines as Oj moves away from this value. 



was calculated (including the full IC cross section) for a jet 
with the parameters of the best fit shown in Fig. [2] but as- 
suming a power law distribution p — 3 from 7 C = 100 up to 
7c, cutoff ~ 3 10 3 . (A p = 3 slope is expected for a steady state 
distribution of electrons injected with the canonical p = 2 
power law in the presence of strong Thompson IC cooling.) 
The >100 MeV lightcurve was indistinguishable from the 
one in Fig. [5J even though the cutoff energy changed signif- 
icantly along the orbit due to Doppler boosting. Hence, the 
results obtained here are likely to extend when more com- 
plex spectral shapes and Klein-Nishina effects are taken into 
account . 



3.3 Jet precession 

The preceding section showed that the jet must be inclined 
in order to obtain good fits to the gamma-ray modulation. 
There is evidence for jet inclination in Cyg X-3 as well as 
other microquasars l|Maccaronel [2002) . An inclined jet is 
likely to undergo precession on a timescale longer than the 
orbital period. There is currently no evidence for or against 
jet precession in Cyg X-3. Here, jet precession will manifest 
itself as a change in the gamma-ray modulation since Oj will 
sample the full range from to 27r in a full precession. Both 
the shape and amplitude are affected as shown in Figure [4] 
The peak flux phase and amplitude can vary dramatically 
from one precession phase to another. 

The Fermi data already show a hint for a change in 
the phasing of the modulation between the two epochs 
during which Cyg X-3 was detected. In addition, the first 
reported detection of Cyg X-3 at 100 MeV from SAS- 
2 showed a gamma-ray orbital modulation correlated (in- 
stead of roughly anti-co rrelated) with the X-ray modu- 
lation (|Lamb et al.l 1 19771) . Later observations by Cos B 
and EGRET fail ed to re-detect the source unambiguously 
|Mori et al.lll997l ). A possible explanation is that the jet ori- 
entation had changed in between these observations. Future 
Fermi observations of Cyg X-3 may find a different mod- 
ulation lightcurve or may actually fail to detect the source 
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because of its low flux, even though Cyg X-3 shows the right 
radio and X-ray state. 

The comparison between gamma-ray lightcurves can 
serve as a very powerful diagnostic of the jet geometry. For 
instance, in the microblazar models discussed in §3.2, the 
near perfect alignment of a jet with the line-of-sight and the 
high P means that the gamma-ray flux is detectable only 
during the very short interval in precession phase where it 
is Doppler boosted. The gamma-ray flux will be deboosted 
most of the time — so that the Fermi and AGILE detections 
would have required very special circumstances. 



4 CONCLUSIONS 

The orbital modulation of the MOO MeV flux from Cyg X-3 
can be very well fitted by a simple-minded model in which 
the emission is due to HE electrons up-scattering stellar pho- 
tons. The HE electrons are situated in two symmetric loca- 
tions in a relativistic jet with an arbitrary orientation. 

The fitting procedure reveals that the jet is necessar- 
ily inclined to the orbital plane normal. The most likely 
value is close to the line-of-sight (</>j ~ i, in agreement 
with the conclusion s base d on radio imaging of the jet 
(Mio duszewski et al.l l200lh . The HE electrons cannot be 
close to the compact object. They are outside of the sys- 
tem at distances of at least one orbital separation, possibly 
up to lOd. IC scattering of accretion disc photons is then ir- 
relevant. If the compact object in Cyg X-3 is a neutron star, 
the required power in HE electrons is a significant fraction 
of the Eddington luminosity. For a black hole, because of the 
lower system inclination implied, the power required can be 
as low as 10~ 5 Z/Edd- These conclusions appear robust even 
when more complex electron distributions and the full IC 
cross-section are taken into account. Precession can be ex- 
pected from an inclined jet. It should cause a change in the 
shape and amplitude of the gamma-ray modulation in the 
future. 

The IC cooling timescale is i ic « 0.5(7 e /10 3 ) _1 (ii/d) 2 
seconds (scaled to the orbital separation d and for orbit Ol). 
The size of the gamma ray emitting region is roughly s w 
/3ct ic , giving s/R < 0.04/3(7e/10 3 ) -1 (i?/d) when scaled to 
R. Hence, the assumption that the emission in the Fermi 
energy range is localised holds up to distances « lOd from 
the star. Cooling slows down at lower energies and electrons 
emit synchrotron radio beyond the 7-ray emission zone on 
much larger scales. 

The 7-ray emission zone could be related to electron ac- 
celeration at a recollimation shock as the jet pushes its way 
through the stellar wind. The jet is initially over-pressured 
compared to its environment. It expands freely until its pres- 
sure pj matches that of the environment p e . Here, p e is the 
ram pressure of the supersonic wind pwV^. The jet pres- 
sure is Pj ~ Lj /(tycQ 2 1 2 ) where Lj is the jet power, Q is 
its opening angle and I is t he distance along the jet (e.g. 
iBednarek fc Protheroelll997fl . The pressures equilibrate at 

i-o-se-^AcM (7) 

with Lj = 10 38 erg s" 1 , M w = 1O _5 M yr -1 and v w = 
1000 km s" 1 . A jet recollimation shock forms beyond I. The 
shock crosses the jet axis after a further distance of order I 



when the external pressure is constant (|Stawarz et alj|2006h . 
This is roughly the case here since the jet does not extend 
very far from the system and the dependence of p w with / 
remain shallow (unless it is pointed directly away from the 
star). The location is consistent with the values of H derived 
above, suggesting this is where jet kinetic or magnetic en- 
ergy is channeled into particle acceleration. This should be 
verified by calculations taking into account the non-radial 
nature of the jet-wind interaction. The shock occurs in the 
wind only because M w is very large (WR star) and the orbit 
very tight. Most microquasar jets will actually break out of 
the immediate vicinity of the system and interact much fur- 
ther away when their pressure matches that of the ISM. Any 
HE particles there will find a much weaker radiation envi- 
ronment and will be less likely to produce a (modulated) IC 
gamma-ray flux detectable by Fermi or AGILE. 

The emerging picture is that of a jet launched around a 
black hole, with a moderate bulk relativistic speed, oriented 
not too far from the line-of-sight, interacting with the WR 
stellar wind to produce a shock at a distance of l-10d from 
the system, where electrons are accelerated to GeV energies 
and upscatter star photons. 
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